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ORIGINAL ARTICLE 

Increased binding of MeCP2 to the GADl and RELN 
promoters may be mediated by an enrichment of 5-hmC in 
autism spectrum disorder (ASD) cerebellum 

A Zhubi\ Y Chen\ E Dong\ EH Cook^ A Guidotti^ and DR Grayson^ 

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by symptoms related to altered social interactions/ 
communication and restricted and repetitive behaviors. In addition to genetic risk, epigenetic mechanisms (which include DNA 
methylation/demethylation) are thought to be important in the etiopathogenesis of ASD. We studied epigenetic mechanisms 
underlying the transcriptional regulation of candidate genes in cerebella of ASD patients, including the binding of IVleCP2 (methyl 
CpG binding protein-2) to the glutamic acid decarboxylase 67 (GADl), glutamic acid decarboxylase 65 (GAD2), and Reelin (RELN) 
promoters and gene bodies. Moreover, we performed methyl DNA immunoprecipitation (MeDIP) and hydroxymethyl DNA 
immunoprecipitation (hMeDIP) to measure total 5-methylcytosine (5-mC) and 5-hydroxymethylcytosine (5-hmC) in the same 
regions of these genes. The enrichment of 5-hmC and decrease in 5-mC at the GADl or RELN promoters detected by 5-hmC and 5- 
mC antibodies was confirmed by Tet-assisted bisulfite (TAB) pyrosequencing. The results showed a marked and significant increase 
in l\/leCP2 binding to the promoter regions of GADl and RELN, but not to the corresponding gene body regions in cerebellar cortex 
of ASD patients. Moreover, we detected a significant increase in TET1 expression and an enrichment in the level of 5-hmC, but not 
5-mC, at the promoters of GADl and RELN in ASD when compared with CON. Moreover, there was increased TET1 binding to these 
promoter regions. These data are consistent with the hypothesis that an increase of 5-hmC (relative to 5-mC) at specific gene 
domains enhances the binding of MeCP2 to 5-hmC and reduces expression of the corresponding target genes in ASD cerebella. 
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INTRODUCTION 

Idiopathic 'Autism Spectrum Disorder' (ASD) is a neurodevelop- 
mental disorder characterized by symptoms associated with social 
interactions, communication skills and restrictive or repetitive 
behaviors. The clinical presentation of autism has been reported 
to overlap with several related neurodevelopmental disorders 
including Rett syndrome (RTT) and methyl CpG binding protein-2 
(MeCP2) duplication syndrome.^'^ Altered MeCP2 levels due to 
gene deletions or mutations are responsible for the majority of 
RTT cases.^ Furthermore, an increase in MeCP2 gene dosage found 
in MeCP2 gene duplication cases'^ impacts a critically important 
epigenetic mechanism, hence altering brain development and 
leading to an autistic-like phenotype. 

MeCP2 is a member of the methyl-binding domain family of 
proteins that is highly expressed in CNS neurons and is particularly 
enriched in GABAergic interneurons.^'^ MeCP2 binds with high 
affinity to modified CpG dinucleotides and is an essential 
epigenetic regulator operative during human brain 
development.^'^ The functional importance of MeCP2 in the 
regulation of synaptic and neuronal plasticity, development of 
motor skills, cognitive and social behavior is well documented.^ 
While idiopathic ASD is rarely associated with MeCP2 mutations or 
duplications, reduced MeCP2 levels are reported in brains of ASD 



patients.^ This suggests that in addition to gene dose effects, 
epigenetic mechanisms may also regulate the expression of 
MeCP2.^"" 

The role of MeCP2 in regulating transcription is more complex 
than previously proposed and it also appears to act as a positive 
modulator of gene expression, rather than strictly as a transcrip- 
tional repressor.^ Whether MeCP2 acts to repress or activate genes 
likely depends on post-translational modifications (e.g. phosphor- 
ylation, acethylation, ubiquitination) and/or cytosine modification 
status (i.e. 5-hydroxymethylcytosine (5-hmC) or 5-methylcytosine 
(5-mC)) of associated target genes.^^'^^ Any mechanistic explana- 
tion for the pleiotropic action of MeCP2 (as either a positive or 
negative modulator of transcription) is likely to involve the 
binding affinity of MeCP2 for 5-hmC.^^ In other words, the 
differential action of MeCP2 at 5-hmC- or 5-mC-enriched DNA 
domains may be associated with complex structural transitions, 
which impact promoter accessibility to ancillary factors that are 
either positive or negative. 

CpG and non-CpG DNA methylation proximal to promoters or 
within gene bodies has been generally regarded as a highly stable 
epigenetic mark that ensures gene expression homeostasis and 
maintains neuronal phenotype integrity.^^ However, recent 
studies suggest that this epigenetic marking is highly dynamic. 
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It Is now believed that steady-state DNA methylation is 
maintained by the opposing action of DNA-methyltransferases 
(DNMT) and an active DNA-demethylation pathway.^ ^ DNA 
demethylation involves the initial hydroxylation of 5-mC to form 
5-hmC by members of the TET protein family.^^'^^ The adult 
mammalian brain expresses high levels of 5-hmC reaching values 
of approximately 40% of methylated CpG sites in Purkinje 
cells.^^'^^ Compared with the other brain areas, the cerebellum 
of RTT patients contains enriched amounts of 5-hmC.^^'^° Recent 
studies show that there is an inverse relationship between 5-hmC 
(enriched) and 5-mC (depleted) abundance in the gene body 
domains of genes expressed in the cerebella of RTT patients.^ 

Because the levels of 5-hmC in relation to 5-mC have not been 
extensively studied in the cerebella of ASD patients, we sought to 
further examine the role of l\/leCP2 in reading epigenetic marks (5- 
hmC/5-mC) of active target genes whose expression is altered. We 
focused on the cerebellum which is anatomically, histologically 
and functionally a well-defined brain structure that has been 
associated with ASD pathogenesis.^^ "^^ One well-established 
finding in postmortem studies of ASD patients is a decrease in 
glutamic acid decarboxylase 67 (GAD1), glutamate decarboxylase 
65 (GAD2), and reelin (RELN) expression in the cerebellum.^^"^^ 
Hence we selected GAD1, GAD2, and RELN as target genes for 
these epigenetic studies. We report that l\/leCP2 binding at the 
GAD1 and RELN promoters is increased and that this increase is 
associated with an enrichment of local 5-hmC levels relative to the 
levels of 5-mC in the cerebella of ASD patients. 



MATERIALS AND METHODS 

Subjects and cerebellar collection 

Blocks of cerebellar cortex from 10 CON and 10 ASD were obtained frozen 
at -80 °C from the Harvard Brain Tissue Resource Center, McLean Hospital 
(Belmont, MA, USA) with approval of Autism Speaks (Autism Tissue 
Program). For each sample, we received a coronal section (0.5 cm block) of 
cerebellar cortex inferior to the horizontal fissure. From the available 
neuropathology reports at the level of where the surface was cut, there 
were no signs of infarction, hemorrhages or inflammatory lesions, with the 
exception of one ASD subject, which had a large early intermediate stage 
infarct. Subjects were roughly matched for age, sex, PMI and brain mass. 
The demographic information associated with each group of subjects is 
listed in Supplementary Table 1. Patients with Asperger's syndrome, 
Fragile-X syndrome, RTT, pervasive developmental disorder not otherwise 
specified, and 15q11-q13 duplication were not included. In addition, all 
control cases were free of neurological disorders, seizures, mental 
retardation, dementia, and so on based on medical records. 

RNA isolation and quantification of mRNA using real-time PCR 
mRNA extraction, purification and measurements of mRNA expression 
were conducted using a previously published protocol.^^ mRNAs 
corresponding to DNMT1, MeCP2, TET1, GAD67, GAD65 and RELN were 
measured. ACTB and GAPDH were used as internal controls for sample 
normalization (see Supplementary Table 2). Each target was run in 
duplicate and the two housekeeping genes were run in parallel. Values 
were calculated as relative abundance to the mean of the two 
housekeeping genes after normalization. 

Western blot analysis 

For protein quantification we conducted measurements as described in 
detail elsewhere.^^ Anti-MeCP2 polyclonal antibodies (1:2000 dilution, 
Diagenode, Denville, NJ, USA) were used to detect MeCP2 protein (-75 
kDa). Changes in MeCP2 protein levels in ASD vs CON were determined by 
expressing the data as a percentage of the CON. 

Chromatin immunoprecipitation assays 

We performed chromatin immunoprecipitation (ChIP) assays based on 
protocols as previously described.^^ The percentage of immunoprecipi- 
tated DNA were calculated using the following: % (gene- IP/total 
input) = 2^^^^^°°/° input) - 3.32) -ct(gene- IP) ^^QQo/o, jETI (Zymo, Irving, CA, 



USA) and DNMT1 antibodies (Abeam, Cambridge, MA, USA) were 
previously shown to bind to single immunoreactive bands on western 
blots.^^'^^ 

Methylated and hydroxymethylated DNA immunoprecipitation 
(MeDIP and hMeDIP) 

Methylated/hydroxymethylated DNA immunoprecipitation experiments 
were performed as previously described^^ using 5-mC and 5-hmC 
monoclonal antibodies (Diagenode, Denville, NJ, USA). The percent 
methylated or hydroxymethylated DNA/ input DNA was calculated as 
described for ChlP. 

Amplicon-specific TAB pyrosequencing 

The distribution of 5-hmC and 5-mC in the GAD1 promoter amplicon was 
determined by a modification of TAB-seq.^° We followed the method 
originally described in detail.^^ Reagents were obtained from Wisegene 
(Chicago, IL, USA). Bisulfite conversion was performed using the EZ DNA 
Lightning Kit (Zymo, Irvine, CA, USA). To determine the efficiency of each 
enzymatic step, spike-in controls were added and independently scored by 
conventional Sanger sequencing. Unmethylated XDNA (cl857 Sam7, 
Promega, Madison, Wl, USA) was methylated in vitro with CpG 
Methyltransferase (New England Biolabs, Ipswich, MA, USA) for a 
methylated control template. To generate a 5-hmC-containing control, a 
275-bp segment corresponding to 1635-1810 bp of pGEMI (Promega) was 
amplified using linearized pGEMI and 5-hydroxymethyl-2'- deoxycytidine- 
5'-triphosphate (Bioline,Taunton, MA, USA ) during PCR amplication.^°'^^ 

DNA from three ASD and three CONs were each immunoprecipitated 
(three IPs/patient) with either 5-hmC or 5-mC antibodies (as described 
above). DNA from each patient sample was separately pooled and 
concentrated. In parallel, three 5-hmC and 5-mC IPs were prepared from 
each patient for bisulfite conversion directly to determine the total number 
of Cs present at each position. Following p-glucosyltransferase protection, 
Tetl -assisted oxidation and bisulfite conversion, DNA was amplified with 
Pfu Turbo Cx DNA polymerase (Agilent, Santa Clara, CA, USA). Percent 5- 
hmC and 5-mC at each position along the GAD1 amplicon was measured 
with pyrosequencing (Epigendx, Hopkinton, MA, USA). The GAD1 bisulfite 
primers (Epigendx, Assay ADS3737) were tested to ensure the absence of 
methylation-based amplification bias. 

Statistical analysis 

The primary test of the difference between CON and ASD for each mRNA, 
protein, ChIP, hMeDIP and MeDIP assays was an independent two sample 
f-test for equal and unequal variance as appropriate. Additional analyses 
were conducted using analyses of variance and covariance controlling for 
factors such as gender, PMI, pH, presence of medication, type of 
medication using PASW v.18 software (SPSS, Armonk, NY, USA). Relation- 
ships between mRNA and protein expression variables and ChIP level 
variables were analyzed with Pearson correlations and scatter plots. Two- 
sided probability levels were used for statistical significance (P<0.05), or 
trends (P<0.1). 



RESULTS 

IVleCP2 binding to GABAergic and glutamatergic gene promoters 
and body regions 

As shown by comparing Figures la and b, IVleCP2 binds more 
efficiently to the gene body regions than to the CpG-rich 
promoter regions of RELN, GAD1 and GAD2. The binding of 
l\/leCP2 to these promoter and gene body regions appears to be 
selective, as l\/leCP2 binding is virtually absent (<0.01%) in the 
promoter region of the housekeeping gene GAPDH (-122 bp to 
-1-42 bp). As shown in Figure la, there is a 1.5-2-fold increased 
binding of IVleCP2 to both the GAD1 (P = 0.04) and RELN (P = 0.03) 
promoters in cerebella of ASD compared with CON. In contrast, 
l\/leCP2 binding to the GAD2 promoter region did not show 
significant changes in ASD vs CON (P = 0.32). Furthermore, IVleCP2 
binding to RELN, GAD1 and GAD2 body regions failed to show 
significant changes in ASD (Figure 1b). 

Increased binding of l\/leCP2 protein to the RELN promoter 
correlated with reduced expression of RELN mRNA in ASD 
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Figure 1. MeCP2 binding to RELN, GAD1 and GAD2 pronnoters and gene bodies and correlations with corresponding nnRNA in the cerebella of 
CON and ASD. (a) MeCp2 binding to the promoter regions of RELN ((-220 to +70), (^P = 0.03)) and GAD1 ((-55 to +121), (*P=0.04)) are 
increased in ASD vs CON, while GAD2 ((-1507 to +1310) (P = 0.32)) is unchanged in ASD vs CON. (b) MeCPl binding to gene body regions of 
RELN ((+562 to +763), (P=0.2)), GAD1 ((+656 to +856), (P=0.6)) and GAD2 ((+1293 to +1447), (P=0.9)) are unchanged in ASD vs CON. (c) 
Increased binding of MeCPl to the RELN pronnoter was associated with reduced RELN nnRNA (Pearson r2 = - 0.66, ''P = 0.04, n = 1 0) in ASD. (d) 
Increased binding of MeCPl to the GAD1 pronnoter shows a trend to statistical significance with decreased GAD1 nnRNA expression (Pearson 
r2 = -0.53; P = 0.08, n = 10) in ASD. ASD, autisnn spectrunn disorder; CON, control; GAD1, glutannic acid decarboxylase 67; GAD2, glutannic acid 
decarboxylase 65; RELN, Reelin. 



(Figure 1c). The increased binding of IVleCP2 to the GAD1 
pronnoter was also negatively correlated with GAD1 nnRNA 
expression (Figure Id), however due to the variability and 
relatively snnall nunnber of subjects, this correlation showed only 
a trend towards statistical significance (P< 0.08). 

We tested whether background variables or other possible 
confounding variables may have influenced the binding of l\/leCP2 
to the target genes in the ASD samples. As shown in 
Supplementary Table 1, average values for the confounding 
variables were similar in CON and ASD. Furthermore, ANCOVA did 
not reveal a statistically significant correlation between MeCPl 
binding to either the GAD1 or RELN promoters with changes in 
age (GAD1 vs MeCP2, Pearson = - 0.31, P = 0.19; RELN vs IVleCP2, 
Pearson r^ = -0.39, P = 0.09), sex (GAD1 vs IVleCP2, Pearson 
r2 = 0.28, P = 0.25; RELN vs IVleCP2, Pearson r^ = -0.14, P = 0.56), 
PMI (GAD1 vs IVleCP2, Pearson r^ = -0.3, P = 0.23; RELN vs IVleCP2, 
Pearson /^ = -0.33, P = 0.17) in CON and ASD patients considered 
either together or separately. In addition we did not find a 



statistically significant correlation between PMI and GAD1 
promoter 5-hmC (Pearson i^ = - 0.22, P = 0.45) or 5-mC enrichment 
(Pearson = - 0.033, P = 0.914), and RELN promoter 5-hmC 
(Pearson /^ = -0.39, P = 0.15), or 5-mC enrichment (Pearson 
r^ = -0.02, P = 0.94) in CON and ASD patients. Race was a 
confounding factor for which we had insufficient information 
and so was not included in the above analysis. This was 
particularly the case in the CON group (see ethnicity. 
Supplementary Table 1). In a subpopulation of ASD subjects with 
seizure comorbidity (A/ = 4), including one subject with signs of 
cerebellar infarct, the levels of mRNAs encoding MeCPl, GAD1, 
RELN; IVleCP2 binding to GAD1 and RELN; and 5-hmC enrichment 
of the GAD1 and RELN promoters, were each similar to the 
average values obtained from the entire ASD population. 

We found no differences in IVleCP2 binding between ASD 
subjects who were treated (antipsychotics, antidepressants or 
mood stabilizers (n = 7)), duration of treatment, and the few ASD 
(n = 3) that apparently were non-treated at the time of death. 
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Figure 2. MeCPl mRNA and protein expression in CON vs ASD. (a) MeCPl mRNA expression is increased in ASD vs CON (^P = 0.01). Values 
were corrected for expression of housekeeping gene GAPDH mRNA. All values are expressed as means ±s.e.m. (b) MeCPl protein levels fail to 
increase in ASD samples relative to CON (P = 0.45, n = 10). MeCPl protein levels in ASD were expressed as percentage relative to CON (100%). 
(c) Representative immunoblots show a major band of 75 kDa for MeCPl and 42 kDa for p-actin. ASD, autism spectrum disorder; CON, control. 
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Figure 3. Connparison of 5-hnnC and 5-nnC enrichnnent in two different gene regions of the RELN, GAD1, GAD2 and correlations between 5- 
hnnc enrichnnent in the pronnoter regions of RELN and GAD1 with the respective nnRNA levels in the cerebella of CON and ASD. (a) Increased 
enrichment of 5-hmC at the promoters of RELN ((-220 to +70 bp) (*P = 0.026)) and GAD1 ((-55 to +121 bp), (*P = 0.047)), but not GAD2 ((-1507 
to -1310 bp), (P = 0.36)) in cerebella of ASD vs CON. No changes in the enrichment of 5-hmC at the gene bodies of RELN ((+562 to +763 bp), 
(P=0.22)), GAD1 ((+562 to +763 bp), (P = 0.53)) and GAD2 ((+1293 to +1447 bp), (P = 0.46)) in cerebella of ASD vs CON. (b) No changes in the 
enrichment of 5-mC at the promoter of RELN ((-220 to +70 bp), (P= 0.91 )), GAD1 ((-55 to +1 21 bp), (P = 0.67)) and GAD2 ((-1 507 to -1 31 0 bp), 
(P = 0.36)) in cerebellum of ASD vs CON. There were no changes in the enrichment of 5-mC at body regions of RELN ((+562 to +763 bp), 
(P = 0.16)), GAD1 ((+656 to +856), (P=0.11)) and GAD2 body ((+1293 to +1447 bp), (P=0.46)) in cerebella of ASD vs CON. (c) Statistically 
significant correlation of 5-hmC in the promoter of GAD1 ((-55 to +121 bp), (Pearson r2 = -0.40, *P = 0.036)) and RELN ((-220 to +70 bp), 
(Pearson r2 = -0.65, ^P = 0.01)) with corresponding mRNA levels in ASD samples. No significant correlations between 5-mC content in the 
promoters of GAD1 ((-55 to +121 bp), (Pearson r2 = -0.14, P = 0.71)) and RELN ((-220 to +70 bp) vs ASD (Pearson r2 = -0.10, P=0.75)) with 
mRNA levels were evident in ASD (vs CON) samples. ASD, autism spectrum disorder; CON, control; GAD1, glutamic acid decarboxylase 67; 
GAD2, glutamic acid decarboxylase 65; RELN, Reelin. 



Hence, the study shows that a diagnosis of ASD is a significant 
predictor for increased l\/leCP2 binding to the GAD1 and RELN 
promoters in cerebellar cortex of ASD patients. 



IVleCP2 mRNA and protein levels in cerebella of ASD patients 
As shown in Figure 2a, the expression of IVleCP2 mRNA was 
increased (by -70%) in ASD cerebella when compared with CON. 
Overall, correlation analysis showed a positive correlation between 
MeCP2 mRNA level and binding of IVleCP2 protein to the RELN 
(Pearson r^ = 0.54, P = 0.04) and GAD1 promoters (Pearson 
r^ = 0.59, P = 0.01). In contrast, there were no correlations between 
IVleCP2 mRNA expression and binding of IVleCP2 protein to the 
RELN, GAD1 and GAD2 gene body regions. 

When measured in cerebellar cortical homogenates, the level of 
IVleCP2 protein was also increased in six and decreased in four of 
ten ASD vs CON. Hence, the increase in mean l\/leCP2 protein 
levels did not reach statistical significance (Figure 2b). 



Enrichment of 5-hmC, but not of 5-mC at RELN and GAD1 
promoters 

As shown in Figure 3a, the levels of 5-hmC detected by 5-hlVleDIP, 
were significantly enriched at the GAD1 and RELN promoters of 
ASD patients whereas 5-mC (Figure 3b) failed to change. 
Interestingly, there is a statistically significant inverse correlation 
between the 5-hmC content at the GAD1 and RELN promoters and 
their corresponding mRNA levels (Figure 3c) in ASD patients. A 
similar correlation was not observed between the GAD1 and RELN 
promoter 5-mC content and the corresponding mRNA levels 
(Figure 3c). Furthermore, there was no significant enrichment in 
the level of 5-hmC and 5-mC in the gene body regions of RELN, 
GAD1 and GAD2 in ASD vs CON (Figures 3a and b). The levels of 5- 
hmC or 5-mC in the gene body regions of GAD2 were very low 
and virtually not distinguishable from background (Figures 3a and 
b). This result indicates a degree of specificity of the observed 
changes that appear restricted to the promoter regions of the 
genes analyzed. 

In order to validate the results obtained from the MeDIP and 
hMeDIP assays, we chose three representative samples for further 
analysis of the GAD1 promoter enrichment by Tet-assisted 
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bisulfite pyrosequencing with base resolution. Our study shows 
(Figure 4) that in CON subjects, 5-hnnC is preferentially enriched at 
CpG position -2 of the GAD1 promoter (sequence relative to the 
TSS: -54 to +69 bp). In the sanne sannples, the levels of 5-nnC is low 
but detectable at positions -2, -1 and 3 and absent in other CpG 
sites studied. In ASD patients (as compared with CON), 5-hmC is 
preferentially enriched at CpG position -2 with smaller amounts at 
position 1 and 2. 5-mC is depleted at this site in the ASD samples 
but present at these positions in the CON samples (see Figure 4). 
These methylation changes result in a significant increase in the 
amounts of 5-hmC relative to 5-mC and the overall 5-hmC/5-mC 
ratio at the GAD1 promoter (-54 to +69 bp) is 1.2 in CON and 5.5 
in ASD patient DNA. 

DNMTI and TET1 mRNA levels 

TET1 mRNA was increased in the cerebella of ASD patients and 
this is consistent with the increased binding of TET1 to the GAD1 
and RELN promoters (Table 1a). This is also consistent with the 
increased 5-hmC levels observed at the GAD1 and RELN 
promoters (Figures 3a and b). DNMTI mRNA was virtually 
unchanged in the cerebella of ASD patients as was the binding 
of DNMT1 to GAD1 and RELN promoters (see Table lb). 



DISCUSSION 

By measuring the binding of l\/leCP2 to the CpG-rich promoter and 
body regions of GAD1 and RELN, we unexpectedly detected a 
marked and significant increase in IVleCP2 binding to the 
promoters but not the corresponding body regions of these 
genes in cerebellar cortices of ASD patients. Importantly, 
demographic data analyses of the patient cohort shows that a 



diagnosis of ASD, and not other confounding variables, is a 
significant predictor of increased l\/leCP2 binding to the GAD1 and 
RELN promoters. In addition, increased binding of IVleCP2 is not 
restricted to the cerebellum. In fact, a similar increase of MeCP2 
binding to the GAD1 and RELN promoters is observed in 
prefrontal cortices of the same patients.^^ 

The increased binding of IVleCP2 to the GAD1 and RELN 
promoters could be due to increased levels of IVleCP2, an 
increased affinity of IVleCP2 (or post-translationally modified 
IVleCP2) to DNA, or an increased presence of methylated or 
hydroxymethylated cytosines at the corresponding promoters (or 
a combination of these). We report that, in our cohort, the levels of 
IVleCP2 mRNA and protein are not decreased, but show a 
tendency to increase in the cerebella of ASD patients. However, 
an increase of l\/leCP2 or a post-translationally modified l\/leCP2 
cannot fully explain why the binding of this protein to the GAD1 
and RELN promoters is increased. In fact, increased IVleCP2 binding 
to these promoter regions is also not accompanied by increased 
binding to the respective gene bodies indicating that the levels of 
IVleCP2 (or of a post-translationally modified l\/leCP2) are not a 
limiting factor in determining the distribution of l\/leCP2 at specific 
gene domains (i.e. promoters vs gene bodies). 

A recent study demonstrated that mouse IVleCP2 binds to both 
5-mC and 5-hmC with high affinity.^^ Using MeDIP and hMeDIP 
with modification-specific antibodies, we found an enrichment of 
5-hmC (but not 5-mC) at the GAD1 and RELN promoters and not 
the promoter of GAD2 in ASD patients. In ASD patients that 
express comparable amounts of l\/leCP2 as controls, our data 
demonstrate that there is more l\/leCP2 binding to these genes 
and that the increased binding of IVleCP2 to the GAD1 and RELN 
promoters is likely related to the increased amounts of 5-hmC 
present. 

Both IVleCP2 knockout, conditional-mutant, and IVleCP2 knock-in 
mice have been widely used to model various aspects of IVleCP2 
loss or gain of function in RTT or ASD patients.^"^"^^ Conditional 
loss of l\/leCP2 in vesicular inhibitory amino acid transporter 
expressing GABAergic mouse neurons results in reduced inhibi- 
tory transmitter quantal size, which is accompanied by an -50% 
decrease in the levels of Gadi and Gad2 mRNA and reduced GABA 
immunoreactivity.^ However, according to recent reviews of 
IVleCP2 function in mutant mice, impaired biological function of 
IVleCP2 occurs as the result of loss of expression, underexpression 
and also overexpression.^'^^ Hence, it is not surprising that ASD 
cerebella are associated with an increased l\/leCP2 binding to 
target gene promoters. 

Genome-wide analysis of modified cytosines (5-hmC and 5-mC) 
across DNA from mouse cerebella shows that 5-hmC is enriched in 
actively expressed genes while 5-mC is depleted at these same 
regions.^ ^ The resultant 5-hmC/5-mC ratio is increased and 
functionally coupled with increased transcriptional activity.^^ Our 
data show that while the ratio of 5-hmC/5-mC has a tendency to 
be higher at both the GADI and RELN promoters and gene bodies 
in ASD cerebella (compared with CON), the increased ratio did not 
correlate with increased mRNA expression. In fact, we found that 
the increased 5-hmC content corresponding to the GADI and 
RELN promoters inversely correlated with gene expression. We 
should note that the gene body and promoter regions reported in 
our study are quite small (-200 bps) compared with those used in 
genome wide studies (>50kb).^^ Our results raise the interesting 
possibility that the observed increase in TET1 expression leads to 
elevated levels of 5-hmC, which facilitates increased binding of 
IVleCP2 perhaps preventing transcription. In addition, we demon- 
strate that DNMTI (methylation writer) is not changed in the 
cerebella of ASD patients whereas TET1 is increased significantly in 
these same samples (see Table 1). Moreover, we show that the 
binding of TET1 to candidate target promoters (GAD1 and RELN) is 
increased in ASD cerebella. We previously reported that TET1 
mRNA is increased in the cortex of SZ patients and that this 
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Table 1. Summary of mRNA levels and DNA binding of TET1 (a) and DNMTI (b) to the RELN and GAD1 promoters in cerebella of ASD vs CON 



(a) Chip (% of DNA input) (b) CHIP (% of DNA input) 



TETl relative mRNA TETl binding to RELN TETl binding to GAD! DNMTI relative mRNA DNMTI binding to RELN DNMTI binding to GADl 



CON 0.021 ± 0.01 0 0.01 1 ± 0.021 0.070 ± 0.01 1 0.76 ± 0.25 0.20 ± 0.01 0 0.1 1 ± 0.020 

ASD 0.052 ±0.01 0 0.031 ± 0.064 0.1 00 ± 0.01 3 0.60 ± 0.1 1 0.1 6 ± 0.030 0.1 3 ± 0.020 

P-value 0.01 0.01 0.046 0.6 0.5 0.6 



Abbreviations: ASD, autism spectrum disorder; ChIP, Chromatin ImmunoPrecipitation; CON, control. % of DNA input was calculated as 2^'^*^^°°''° input)- 3.32) -ct 
(gene- IP) ^ 100% where IP is the amount of specific gene-immunoprecipitated DNA and input is the amount of the same gene in the starting material. Relative 
mRNA levels were calculated relative to the amounts of GAPDH mRNA. 



correlates with increased levels of 5-hnnC at the GADl 
pronnoter.^^'^^ 

Connparisons of wild type and l\/leCP2 null (knockout) nnice 
showed that IVleCP2 has no role in establishing the levels of 5-hnnC 
and that regional 5-hnnC/5-nnC ratios occur by a nnechanisnn 
independent of l\/leCP2 binding or expression. Instead, the 
observed increase in TETl expression in ASD cerebella is 
associated with an enrichnnent of 5-hnnC at selected pronnoters, 
including GADl and RELN. Hence, our data suggest that the 
increased binding of l\/leCP2 to the GADl and RELN regulatory 
donnains nnay be nnediated by an enrichnnent of 5-hnnC. Taken 
together, these findings extend a role for IVleCP2 from IVleCP2 
duplication syndronne to idiopathic ASD patients. Moreover, these 
data underscore the innportance of regional variations in 5-hnnC 
and 5-nnC content in the regulation of GADl and RELN expression. 
It seenns plausible to suggest that these regional variations nnay 
lead to alterations in the dissociation of MeCP2 fronn specific gene 
targets in ASD cerebella. Collectively, we believe that recent 
studies^^'^° (including our own) warrant additional investigations 
into the role of 5-hnnC in nnodulating differential transcription in 
ASD brains. 
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